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ABSTRACT 

(N 

Due to the co-evolution of supermassive black holes and their host galax- 
ies, understanding the mechanisms that trigger active galactic nuclei (AGN) are 
imperative to understanding galaxy evolution and the formation of massive galax- 
ies. It is observationally difficult to determine the trigger of a given AGN due 
to the difference between the AGN lifetime and triggering timescales. Here, we 
utilize AGN population synthesis modeling to determine the importance of dif- 
ferent AGN triggering mechanisms. An AGN population model is computed by 
combining an observationally motivated AGN triggering rate and a theoretical 
AGN light curve. The free parameters of the AGN light curve are constrained by 

^ minimizing a \ 2 test with respect to the observed AGN hard X-ray luminosity 

function. The observed black hole space density, AGN number counts, and X- 

iy) ! ray background spectrum are also considered as observational constraints. It is 

found that major mergers are not able to account for the entire AGN population. 
Therefore, non-merger processes, such as secular mechanisms, must also trigger 
AGN. Indeed, non-merger processes are the dominant AGN triggering mecha- 
nism at z < 1—1.5. Furthermore, the shape and evolution of the black hole mass 

^ function of AGN triggered by major mergers is intrinsically different from the 

shape and evolution of the black hole mass function of AGN triggered by secular 
processes. 

Subject headings: galaxies: active - - galaxies: Seyfert - - galaxies: quasars: 
general — X-rays: diffuse background 



1. Introduction 

In less than a decade it was determined that not only do all massiv e galaxies har- 



bor a supermassive black hole at their center (JKormendy fc Richstondll995l ). but also that 
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the e volution of the galaxy is intrinsically linked to the growth of the central black hole 



[e.g., iMagorrian et al.l Il998l ; iTremaine et al.l 120021 ) . However, the mechanism responsible 
for the co-evolution of a supermassive black hole and its host galaxy is still under inves- 
tigation. As accreting supermassive black holes, generally referred to as active galactic 
nuclei (AGN), can, over their lifetime, r adiate an amount of energy comparable to the 



1999: Alexander et al. 



bindi ng energy of their host galaxy (e.g., ISilk fc Reed Il998l ; iFabian 

20101 ) , it is believed that this AGN feedback will affect star formation in the host galaxy 



ag. 



Fabian 



2010 



1999; 



Hopkins et al. 



lHambrick et al 



2011 



2006a l iLagos et al. 



20081 : iTrichas et al.ll200a Bertone et al. 



Hocuk fc Spaansl l201ll ). Galaxy wide processes can also af- 



fect the central supermassive black hole. For example, theoretical and observational evi- 
denc e suggests that AGN can be triggered by major mergers of massive gas rich galaxies 



Sanders et al.lll988l : iHernquist et al.lll98a ICarlberdll99d : Kauffmann fc Haehneltll2000 



Hopkins et al.l l2006aj ) or by secular processes not connected to major mergers, such as su- 
pernova winds, stellar bars, cold-flow a ccretion, interactions with other massive galaxies or 



satellite galaxies, or minor mergers (e.g . .ICrenshaw et al. 



Vittorini et al.l2005: Davies et al. 2007: 



2010l : iBournaud et al.ll201ll ; lOrban de Xivry et al. 



Pierce et a 



2007 



20111 ). 



2004 



2003; Korme ndv fc Kennicutt 
Chen et al.ll2009t iKumar fc Johnson 



Simulations show that mergers of gas rich galaxies cause g as and dust to lose angula r 



momentum and fall into the central regions of the galaxy (e.g., iBarnes fc Hernquistlll996l ). 
The resulting nuclear gas reservoir will be consumed by star formation in the hos t bulge 
and accretion flows onto the resident supermassive black hole (IHopkins et al.ll2006al Jbl). Ac- 
cording to the galaxy evolution model explored by iHickox et al.l (120091 ). AGN feedback will 
then evolve the AGN host galaxy across the color- magnitude diagram (CMD) from the mas- 
sive end of the blue cloud, through the green valley, and onto the red sequence (but see 
Schawinski et al. 2009 and Cardamone et al. 2010). Secular processes, such as supernova 
explosions, stellar bars, minor mergers and interactions, will also release angular momentum 
from gas and dust leading to a similar reservoir of gas in the central region of the galaxy (e.g., 
Crenshaw etalll2003l : Bavies et al.ll2007l : Ichen et al.ll2009[ JOrban de Xivrv et al.llioilh . De- 
spite mergers and secular mechanisms both leading to the accumulation of gas and dust 
deep in the galactic potential well, mergers and secular processes have very different galaxy 
wide effects and roles in galaxy evolution. Major mergers are violent processes which, as 
with the Antennae galaxies, can destroy g alactic structure, causing th e merger remnant to 
relax into a bulge-dominated system (e.g., IHopkins fc Hernquistl 120061 ) . Secular evolution, 



how ever, is not likely to dist urb the morphology or large-scale structure of the host galaxy 



sec 



Georgakakis et al.ll2009l ). Despite the systemic differences between galaxy evolution and 



AGN activity due to major mergers and secular processes, it is difficult to observationally 
determine which mechanism is responsible for a given AGN; tidal tails are often faint, bulge- 
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dominated systems can be re-triggered by secular processes, and the timescales for AGN 
activit y tend to be longer tha n the timescales for a merger remnant galaxy to dynamically 
relax (ISchawinski et al.ll2010al ). 



Despite the observational difficulties in determining how a particular AGN has been 
triggered, there is observational evidence that major mergers are not the dominant AGN 
triggering mechanis m, at least not a t z < 2. By determining the Sersic indices of massive 
galaxies at z — 2-3, IWeinzirl et al.l (120111 ) found that ~65% of AGN hosts at this redshift 
range have Sersic indices indicati ve of disky morpholog ies. Observations also show that 
in th e redshift ranges z = 1.5-3 ( ISchawinski et al.ll201ll ) and z = 1.5-2.5 ( jKocevski et al. 



20121 ). the majority o f moderate luminosity AGN are hosted by disk galaxies. By considering 
galaxy pairs at z < 2, IWilliams et al.l ( 1201 ll ) find that the wet major merger rate is too low to 
account for the majori ty of AGN activity at z < 2. When investigating AGN in zCOSMOS, 
Silverman et al.l ( 1201 ll ) found that ~20% of moderate luminosity AGN at 0.25 < z < 1.05 



are in close pairs, and thus were likely triggered by galaxy interactions. ICisternas et al. 



( 1201 ll ) argue that AGN are as likely as quiescent galaxies to show signs of a recent merger 
at z < 1, and thus mer gers cannot be connected to AGN activity at this redshift range. 
Georgakakis et al.l (120091 ) investigate the morphology of AGN hosts and find that AGN hosted 
by disk galaxies contribute ~25% of the AGN luminosity density at z « 0.8. They conclude, 
therefore, that a large fracti o n of z ~ 1 AGN are triggered by mechanisms unrelated to major 



mergers. 



Cardamone et al.l (120101 ) found that when AGN host galaxy colors are corrected 
for dust extinction, there is a bi- modality of AGN host galaxy colors at z ~ 1, suggesting 
that there are two modes of AGN activity at this redshift range. W hen m o delin g the AGN 
population by considering mergers of massive dark matter halos, IShenl (120091 ) find that 
secular AGN acti vity is necessary at z < 0.5 to account for the observed AGN population. 
Koss et all (boioh recently found that ~20% of Swift/BAT AGN at z < 0.05 are hosted by 
galaxies with disturbed morphologies indicative of a recent major merger and an additional 
6% of Swift/BAT AGN are in close pairs suggesting these AGN were triggered by galaxy 
interactions. Thus, it appears that both major mergers and secular processes must contribute 
to the AGN activity observed at z < 2. If secular triggers dominate the AGN population at 
z < 2, current understanding of the stochastic fueling of secul ar mechanisms is insufficien t 
to explai n the high luminosity of quasars observed at z < 2 (IHopkins fc Hernquistl 120061 ) . 
However, iBournaud et al.l (120 111 ) recently showed that, at least at high redshift, it is possible 
for cold-flow accretion to trigger moderate luminosity AGN with occasional bright episodes. 
Thus, phenomenological models of the AGN population may provide an important tool for 
understanding the importance of merger and secularly triggered AGN at various redshift 
and luminosity ranges. 



Here, a model of the AGN population is calculated using an empirically motivated space 
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density of AGN triggered at each redshift and a theoretical AGN light curve. We consider 
whether major mergers or secular mechanisms alone can account for the observed AGN hard 
X-ray luminosity function (HXLF), AGN number counts, black hole space density, and the 
X-ray background (XRB) spectrum. In Section [2] the details of the AGN model, including 
the triggering rate, light curve, and black hole mass functions utilized, are described. Section 
[3] explains the calculations completed to compare the model against observations. We then 
consider the model results if AGN are a single population (Section HJ) or two populations 
(Section [5J). The results are discussed and summarized in Sections E] and [7J A ACDM 
cosmology is assumed with Hq = 70 km s _1 Mpc -1 and Q\ = 1.0 - Qm = 0.7. 



2. The AGN Population Model 

In order to model the evolution of the AGN HXLF, three ingredients are necessary. To 
determine the space density of triggered AGN, an empirically based major merger rate is 
used. Once triggered, the AGN Eddington ratio, A = L bo i/L Edd , where L bo[ is the bolometric 
luminosity and Lem is the Eddington luminosity, is evolved using a t heoretically mo t ivated 



light curve. The active black hole mass function (ABHMF) and the iMarconi et al.l (120041 ) 
bolometric correction are then used to convert from Eddington ratio to 2-10 keV luminosity, 
Lx- These three ingredients are described in detail below. 



2.1. Triggering Rate 

In order to determine the space density of AGN triggered at each redshift, the space 
density of gas rich massive galaxies at redshift z must be calculated. First, the minimum 
stellar mass of a massive galaxy, M™ in , at z is derived by parametrizi ng the median mass o f 



ultra- luminous infrared galaxies (ULIRGS) as a function of redshift ( ITreister et al.ll2010al ). 
which gives 

M™ in (z) = 5 x 10 n (1.0 + z)- L5 Af . (1) 

The space density of massive galaxies at z, N ga i(M* > M™ n (z)), in Mpc -3 , can then be 
calculated by integrating the stellar mass function from M™ m (z) to M™ ax = 10 12 ' 5 M , such 
that 

N gal (M*>M™(z))= a ^f M ;; Z) d\ogM^ (2) 



M™ m (z) 



d\ogM* 



where d<& ga i/d\ogM* is the stellar mass function (SMF) of iPerez- Gonzalez et al.l (120081 ) for 
z < 4. The dependence of the results on the SMF is considered in Sections 14.11 and 14.21 To 
determine the space density of gas rich massive galaxies, the fraction of gas rich galaxies at 
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z, fg jz ), must be determined . Considering observations of the GOODS fields (jDahlen et al. 
20071 ). iTreister et al.l (j2010al ) find that f g (z) can be parametrized as 



/,(*) 



0.11(1 + z) 2 - 
1 



z< 2 
z> 2 



(3) 



Thus the space density, in Mpc 3 , of potential AGN host galaxies at redshift z is f g (z)N ga i(M* > 
M™ n (z)). 

The space density of AGN triggered by a merger at redshift z is then calculated by 
multiplying the space density of potential AGN host gal axies at z by the fracti on of massive 
galaxies which will undergo a merger at z. Following iHopkins et al.l (j2010al ). who derive 
the major merger rate per galaxy per Gyr, d 2 ^/dt dN, from simulations and observational 
constraints, we parametrize d 2 ^/dt dN as 



d 2 ^ 

dtdN 



A(M™ n )(1.0 + z 



,/3(M™ n ) 



where 



and 



A(M* 



0.02 



M, 



2 x 1O 1O M 



1/2' 



Gyr- 1 



P(M m ) = 1.65 -0.15 log 



M ± 



(4) 



(5) 



(6) 



2 x 1O 1O M 0/ 

Here the major merger is said to occur upon the coalescence of the two similarly massive 
galaxies ( IHopkins et al.ll2010bl ). By assuming that every major merger leads to an AGN 
event, with a negligible time delay between the merger and the triggering of the AGN 
activity, the space density of AGN triggered by major mergers, dN merg , at redshift z is 

d 2 V 



dN (z) 

Ujl v mergy^ I 



dNdt 



N gal (M* > Mr n (z)) f g (z) dtMpc' 



(7) 



The rate at which AGN are triggered through secular processes, such as galaxy inter- 
actions, cold gas accretion, and internal disk instabilities, is calculated in a similar manner. 
For AGN triggered by secular proc esses, M™ m (z) = 5 x 10 9 M , in agreement with the 
findings of ISchawinski et al.l (|2010bJ). The fractional rate of massive gas rich galaxies which 



are tr iggered through secul ar processes every Gyr, f sec , is assumed to be constant with red- 



shift. lYamada et al.l (120091 ) found that the fraction of galaxies with M* > 10 11 M Q hosting 
AGN, which are triggered through both secular processes and mergers, is ~0.3. As the frac- 
tion of AGN triggered by secular processes is poorly constrained observationally, and we are 
considering galaxies with M* > 5 x 10 9 M , the constraint f sec ^C 0.3 Gyr -1 is used. The 
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specific value of f sec is set by calculating the predicted AGN HXLF and minimizing a x 2 test 
which compares against the observed HXLF. The space density of AGN triggered by secular 
mechanisms, dN qpr , at redshift z is then 



dN s , 



/, 



„, iC N gal (M* > M™ n {z))f 9 {z)dtM. V c 



(8) 



2.2. AGN Light Curve 

Once the AGN has been triggered, its Eddington ratio is used to parametrize the accre- 
tion as a function of time since t he AGN was tr i ggered. While there is evidence that AGN 
are an intermittent phenomenon, iTreister et al.l (J2010al ) show that quasars can grow most 
of their black hole mass in a single, merger triggered event. Therefore, a single-peaked light 
curve is assumed. 



Based on hydrodynamical simulations, iHopkins fc Hernquistl (120091 ) suggest 



Kt) = [i + (I*IAq) 1/2 



-2//3 



(9) 



where t = t on — £q, where t on is the time since the AGN was triggered, £q = tor} 13 / (2(3 In 10), 
and t , rj, and /3 are fitting parameters which describe the quasar lifetime, maximum Ed- 
dington ratio, and light curve slope, respectiveljo 



2.3. Active Black Hole Mass Function and Its Evolution 



By combining the triggering rate and light curve, the space density of AGN with Edding- 
ton ratio A can be calculated at any redshift z < 4. However, to compute model predictions 
that can be compared to observational constraints, the black hole mass must be used to 
convert Eddin gton ratios into bolometric luminosities. Once the bolometric luminosity is 
computed, the iMarconi et al.l ( 12004] ) bolometric correction is used to determine Lx- The 
distribution of black hole masses is determined by the fractional active black hole mass func- 
tion (ABHMF) at z, which describes the fraction of active black holes at redshift z with 
black hole mass M,. 



■"■For completeness, a wide variety of light curves were tested, including a constant A, a monotonic lin- 
early increasing or decreasing A, a monotonic exponentially increasing or decreasing A, a linearly increasing 
followed by a linearly decreasing A, and an exponentially increasing followed by an exponentially decreasing 
A. However, these light curve models were unable to match the shape of the AGN HXLF. 



-7- 



Two ABHMFs are considered. The first ABHMF is a Gaussian fit to the combined 
type 1 and type 2 AGN ABHMF at z ~ 0.15 observed by iNetzerl (J2009h . The second 
ABHMF considered is a Schechter function with the same slope and critical mass, M crit , 



as the black hole mass function described by iMerloni k Heinzl (J2008h . Both ABHMFs are 
considered over the range log(M™ m /M ) = 5.95 and \og(M™ ax / M Q ) = 10.55. A black hole 
with M. < M™ n would need to accrete at Eddington ratio A > 0.1 to achieve \ogLx > 
41.5. According to the light curve used here, the vast majority of the AGN lifetime is spent 
at A < 0.1. Therefore, black holes with M. < M™ m are expected to only make a very sma ll 
contribution to the observed AGN population ( IMarconi et al.l 12004 ; IMerloni k Heinzl 120081 ). 
Black holes with M. > M™ a x , corr espon d to <C 0.1% of all active black holes according 
to both the IMerloni k Heinzl J2OO8I ) and INetzerl (12009^ ABHMF s, in agreement with the 
maximum black hole mass derived by iNatarajan k Treisterl ( 120091 ) . 



Both the IMerloni k Heinzl ( 120081 ) and INetzerl (120091 ) ABHMFs are derived using scaling 
relationships to det ermine black hole masses. These scaling relationships have an intrinsic 
scatter of ~0.3 dex JPeterson k Bentzlbood : IMerloni k HeinJ200a : lNetzeJl2009l ). This leads 
to an uncertainty in the calculated Lx of a factor ~ 2 x A. For the majority of sources, A 
< 0.1, thus the uncertainty introduced by the ABHMF is a smaller part of the error budget 
than the uncertainty in the major merger rate, which is a factor ~2. 

For both ABHMFs two r edshift evolutions are investigated. The first ABH MF evolution 
uses the continuity equation ( ISmall k Blandfordlll992l ; IMerloni k Heinzl 120081 ) 



dn M {M.,t) d[n M (M.,t)(M(M.,t))} 



Of 



dM 



0. 



(10) 



where um is the ABHMF and (M(M.,£)) is the average accretion rate of black holes with 
mass M. at time t. By integrating this conservation equation forward and backward in time, 
using the observed ABHMF as the boundary condition, the ABHMF can be evolved to any 
redshift z. This evolution assumes that black hole growth occurs through ac cretion and that 



binar y mergers are not the primary mechanism of black hole growth (e.g., IVolonteri et al. 



2003). The second ABHMF evolution considered is based on the observations of Labita et al. 



(120091 ) who found that the maximum black hole mass of the quasar population increases with 

= M^"'(0)(1.0 + ,z) 1,64 . In order to compare t he two considered 



redshift. Thus, Ml 



criti 



evolutions, a power law is fit to M^ nt (z) of the IMerloni k Heinzl (120081 ) ABHMF evolved 



with the continuity equation and it is found that Ml 



crit I 



M. c "*(0)(1.0 + z) a5 . For both 



ABHMFs, the ABHMF is re-normalized at each redshift so that integrating over all black 
hole masses gives 1. 



With the observationally determined AGN triggering rate, theoretical Eddington ratio 



evolution, and the ABHMF in place, the AGN population at z < 4 can be fully modeleco 
To determine the light curve fitting parameters, to, (3, and rj, the AGN population model 
must be compared to observational constraints. 



3. Calculations and Observational Constraints 

In order to constrain the light curve parameters, the AGN population model is compared 
against the HXLF at five different redshifts, the black hole mass density as a function of 
redshift, the Swift/ 'BAT 15-55 keV AGN counts, the 2-10 keV AGN counts, and the XRB 
spectrum. The methods used to calculate these quantities are described below. 



3.1. AGN HXLF 



The HXLF, #j(Lx, z)/d\ogLx, is com puted at five differ en t redshifts, z = 2. 3, 1.2 , 
0.6, 0.3, and 1 Th e HXLFs ob s erved by JUeda et all d2003h. LaFranca et all (120051 ). 
Silverman et al.l ( 120081 ). lAird et al.l (120101 ) . and lUeda et al.l (120 111 ) are each presented in 
different re dshift bins. Theref ore, we compute the predicted HXLF at the central redshift of 
each of the lUeda et al.l (120031 ) redshift bins. 



The bolometric AGN luminosity function at redshift z, d$(Lt, i, z)/d log L, is computed 
by integrating over the space density of AGN triggered at redshift z t and, at redshift z, have 
black hole mass M, and Eddington ratio A such that the bolometric luminosity is L^. Thus, 

,\ j = / * n M (M.,z)dhgM. / * x (M.,z)dzt, (11) 

UlOgL J M rnin J z 

where log(M™ n /M ) = 5.95, \og(M™ ax / M Q ) = 10.55, z t is the triggering redshift, and 



*x(M.,z) 



§£± X(z,z t )L Edd (M.) = L bol 







otherwise 



(12) 



where dN is either dN merg or dN sec . The iMarconi et al.l ( 12004J ) luminosity dependent bolo- 



metric correction is then used to convert Lh n i to L 



^x- 



As the observed HXLF only includes Compton thin AGN, the Compton thick AGN 
are removed from the calculated HXLF by assuming that one-third of all obscured AGN 



2 AGN triggered at z > 4 make only a minor contribution to the z < 2 AGN population. Therefore, if 
the AGN triggered at z > 4 are included, the results of this study do not change. 
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are Compton thic k in agreement with the fraction of Compton thick AGN necessary for the 
Ueda et al.l (120031) observed HXLF to b e in agreement with th e peak of the XRB at ~ 30 
keV (jDraper fc Ballantynd 120091 . |201Q| ; iBallantyne et al.ll201ll ). The fraction of Compton 
thin obscured sources, / 2 , is assumed t o be a function of luminosity and redshift suc h that 
f 2 oc (1 + z) a (\og L x )~ b , where a = 0.4 feallantvne et al.lbood : iTreister fc Urrvll2006r i and b 
= 4.7. The normalization factor is determined by assuming the type 2 to type 1 AGN ratio 
is 4:1 at z = and \ogL x = 41.5. Thus, the space density of Compton thick AGN also 
depends on L x and, in Mpc~ 3 dex -1 , is (fz/2) d§x/dlogL X - 



According to both the iMarconi et al.l (120041 ) and IVasudevan et al.l (120091 ) bolometric 
corrections, to achieve Lx = 10 46 ' 5 erg s _1 it is necessary to have a black hole wit h mass 
»il0 10 ' 7 - 10 10 ' 9 Mq accreting at itsEddington rate. Both the Gaussian fit to the iNetzer 



(120091 ) and the Merloni fc Heinzl (120081 ) ABHMF, predict the fraction of active black holes 
with mass ~10 11 M m is ~ 0, in agreement with the black hole mass upper limit derived 



by lNatarajan fc Treisterl (120091 ) . This strongly suggests that AGN with Lx > 10 46 erg s 



are accreting at super-Eddington rates. The light curve model used allows for < A < 
1. Therefore, super-Eddington accretion cannot be taken into account. Thus, the models 
presented here will necessarily under-predict the L x > 10 46 erg s _1 observed HXLF data 



points. For completeness, the HXLF data points at L x > 10 
account when performing the x 2 fitting! 



•16 



erg s 1 are still taken into 



3.2. Black Hole Mass Density 



To calculate the black hole mass density, the ISoltanl (119821 ) argument is used. Thus, the 
black hole mass density at redshift z, p»(z), is calculated as 



P*W 



dt 
dz 



dz 



1-e #M,, _ 
-Lboi—r, — ^rflogL 



ec 



dlogL 



(13) 



where e = 0.1 is the radia tive efficiency, c is th e speed of light, and d$(L,z) /dlogL is t he 
AGN luminosity function (lYu fc Tremainelbood : llVlarconi et aDl2004J : ITreister et aJboilh 0. 

The contribution of Compton thick AGN to the black hole mass density is included. The 



3 Exclusion of the three data points at Lx > 10 46 erg s 1 changes the reduced x 2 values by < 0.1 and 
does not affect the determination of the light curve parameters. 

It is expected that at Eddington ratios A < 10~ 2 , AGN acc retion flows becom e radiatively ineffi c ient a nd 



the radiative efficiency decreases with A (e.g. JCao fc Xull2007l ). If Equation 1 of iMerloni fc Heinj (J2008I ) is 
used to characterize e, the calculated local black hole mass density changes by < 5% since the majority of 
black hole growth occurs during radiatively efficient accretion. 
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pr edicted black hole ma ss density is compared to the local black hole mass de nsity observed 
by IShankar et al.l (120091 ) and the z ~ 2 black hole mass density observed by iTreister et al. 
r |2010bh . 



3.3. X-ray Background Spectrum 



The XRB spectrum model closely follows that described by lDraper fc Ballantynd ( 120091 ) . 
Instead of inputting the observed luminosity function, t he HXLF calcula ted as described in 
Section I3TTI is used. Also, a t orus reflection comp onent (IGilli et al.l 120071 ) is included which 
is computed using "reflion" (IRoss fc Fabian! 120051 ) . The Type 2 fraction, f 2 , is determined 
as in Section I3TT1 The Compton thick fraction, fcr, is defined as the ratio of the number of 
Compton thick AGN to the number of Co mpton thin t ype 2 AGN, and is set to for = 0.5, in 
agreement with the fcr necessary for the lUeda et al.l (120031 ) HXLF to be in agreement with 
the XRB spectrum and the local C ompton thick AGN space density (jDraper fc Ballantyne 
20091 . |2010| ; iBallantyne et al.ll201ll ). The unabsorbed type 1 sources are evenly distributed 
over column densities log Nh = 20.0, 20.5, 21.0, and 21.5. The Compton thin type 2 sources 
are distributed equally over log Nh = 22.0, 22.5, 23.0, and 23.5. The contribution of Compton 
thick AGN to the XRB is included. To do so, it is assumed that Compton thick sources 
evolve similarly to less obscured AGN and are evenly distributed over logNu = 24.0, 24.5, 
and 25.0. 

The AGN number counts in the 2-10 keV and 15-55 keV bands are also calculated. This 
is done by using the same AGN spectra, f 2 , and fcr as in the XRB model described above, 
including the contribution of Compton thick AGN, and by using the HXLF calculated as in 
Section 13.11 



3.4. Summary of Free Parameters 

The light curve fitting parameters — to, which is related to the AGN lifetime, 77, which is 
related to the peak Eddington ratio, and /3, which determines the slope of the light curve — 
are determined by comparing the resulting models against the observed HXLF, evolving 
black hole space density, XRB spectrum, and AGN number counts. The fractional rate of 
massive galaxies which are triggered by secular processes each Gyr, f sec , is also determined 
by comparing against observational constraints. 

The focus of this study is not to fit the AGN light curve parameters, but to investigate 
the contribution to the AGN population of AGN triggered by major mergers and AGN 
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triggered by secular processes. Therefore, a minimization algorithm is not used. Instead, 
the best fit parameters for the AGN light curve are determined by considering models with 
t n = 1.0 x 10 7 , 5.0 x 10 7 1.0 x 10 8 , 2.5 x 10 8 , 5.0 x 10 8 , 7.5 x 10 8 , and 1.0 x 10 9 yrs. 
Hopkins &: Hernquistl ( 120061 ) found that their model results were not very dependent on 77, 
so the values 77 = 0.2, 0.4, 1.0, 2.5, and 3.0 are considered. Steps of 0.05 in the range from 
0.05 to 1.0 are used to determine 0. For f sec the values 0.0, 0.005, 0.01, 0.02, and 0.05 Gyr" 1 
are investigated. The best fit parameters for each model are determined by minimizing a x 2 
test which t akes into account the 127 HXLF data poin t s for z = 2.3, 1.2, 0.6, .3, and 0.1 



presen ted b vlUeda et a 



(12010( 1. and lUedaet al. 



(20031). lLa Franca et al.l (l2005f ). ISilverman et all (120081 ) . lAird et al. 



(J201ll ). By comparing against the observed AGN HXLF, black hole 



mass density, XRB spectrum, and AGN number counts, the model is fully constrained, 
allowing for some conclusions to be reached about the AGN light curve parameters. The 
model sensitivity to the light curve parameters is discussed in Section 16.11 



4. One Population 



4.1. Major Merger Triggered Quasars 



First, it is assumed that all AGN are tri ggered by major merger s, thus, f sec = 0.0. 
The best fit to the HXLF is found using the JMerloni fc Heinz! (bood ) ABHMF with the 



Labita et al. 



(2009) evolution, t = 2.5 x 10 8 yrs, (3 = 0.7, and 77 = 2.5, similar t o tha t 



found bylCaol (120101) . When compared t o the 1 2 7 data points from the lUeda et al.l (12003) 



La Franca et al.l (120051 ). ISilverman et all (120081 ) . lAird et al.l (120101 ). and 



Uedaetal 



(120111 ) 

observed HXLFs in the 5 redshift ranges considered, this model has a reduced x 2 , Xred = 2-4. 
If to is changed by 0.5 x 10 8 yrs, (3 is changed by 0.05, or 77 is changed by 0.5, the resulting 
Xred wm increase by ~ 0.1 - 0.2. The light curve fitting parameters are summarized in Table 
[TJ As shown by the dot-dashed red lines in Figure [TJ this model provides a relatively good 
fit to the observed HXLF at z > 1; however, at z < 1, this model has xled = 2-6 for 85 data 
points. The space density of major mergers at z < 1 is too low to explain the space density 
of observed AGN in this redshift range. Thi s model is in agreement with the observed local 
black hole mass densit y (jShankar et al.ll2009l ) and the z ~ 2 observed black hole mass density 
(ITreister et al.ll2010bl ). However, this model significantly under-predicts the XRB spectrum 
as well as the 2-10 keV and 15-55 keV AGN number counts. Thus, despite major merger 
triggered AGN being able to account for the integrated black hole growth, merger triggered 
AGN cannot account for the space density of the entire AGN population. Specifically, major 
mergers cannot account for the z < 1 AGN population. 



The finding that major mergers are not capable of accounting for the z < 1 AGN 
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population is robust against several assumptions. To test the dependen ce of the merger 
trigge red AGN only model on the SMF, t he z = 1-3- 3 SMF presented by iMarchesini et al. 
(120091 ) and the z — 0.1 SMF pre sented bvlCole et al.l (120011) are used to define a n evolving 



SMF. When using the combined ICole et al.l ( 120011 ) and IMarchesini et al.l ( 120091 ) SMF, the 
merger triggered AGN only model is not able to account for the space density of AGN at z < 
1, with minimum x 2 re d = 3-6. Similarly, if the minimum mass of a potential AGN host galaxy 
(M™ n ) is reduced, the merger only model cannot supply a decent fit to the observational 
constraints. 



4.2. Secularly Triggered AGN 

Next, we consider the scenario where all AGN are trigge red by secular processes. For 
this model, the best fit to the HXLF is found when using the iNetzerl ( 120091 ) ABHMF with 
the continuity equation evolution, t = 2.5 x 10 8 yrs, (3 = 0.8, r\ = 0.4, and f sec = 0.02 
Gyr -1 . This model has Xred = 2-1 for 127 data points. If to is changed by ~ 0.5 x 10 8 
yrs,/3 is changed by 0.05, or f sec is changed by 0.01 Gyr -1 , the resulting Xred will change 
by ~ 0.2. If i] is changed by ~ 0.5 the resulting x 2 red will increase by ~ 0.1. As shown 
by the blue dashed lines in Figure [TJ, this model under-predicts the space density of low 
Lx AGN at z < 1 and under-predicts high Lx sources at z > 1. Furthermore, this model 
significantly under-predicts the local black hole mass density, the XRB, and the 2-10 keV 
number counts. This model slightly ove r-predicts the Sw ift/BAT 15-55 keV number count 
and is inconsistent with the findings of iKoss et all (l201of ) that ~20% of Swift/BAT AGN 
host galaxies have disturbed morphology indicative of a recent major merger. Thus, secular 
processes alone are not sufficient to account for the entire AGN population. Specifically, 
secular processes are not able to account for the AGN population at z > 0.5. When the 
minimum mass o f a potenti a l AGN ho st galaxy is reduced or the evolving SMF defined 
by combining the ICole et al.l (120011 ) and IMarchesini et al.l (120091 ) SMFs is used, the secular 
evolution only model still cannot account for the z > 0.5 AGN population, with minimum 
Xred = 2.0. It is found that AGN triggered by secular mechanisms alone cannot explain the 
observed AGN population. 

Neither merger triggered AGN nor secularly triggered AGN alone are able to account 
for the entire AGN popu lation. Red u cing t he m inimum mass of pote ntial AGN host galaxies 
and using the combined I Cole et al.l (120011 ) and iMoretti et al.l (120091 ) SMF does not change 
this result. The space density of major mergers at z < 1 is too low to account for the 
low redshift AGN population but can account for the local space density of black holes. In 
contrast, AGN triggered by secular mechanisms are not able to account for the build up of 
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the black hole mass density over cosmic time. Thus, we consider the possibility that AGN 
are triggered by both major mergers and secular processes. 



5. Two Populations 

As neither the major merger triggered AGN nor the secularly triggered AGN can alone 
account for the entire AGN population, we now consider a two population model in which 
AGN are assumed to be triggered by both secular processes and major mergers. First, the 
major merger triggered AGN light curve is determined by finding the lowest x 2 fit to the z 
= 2.3 HXLF using only the major merger triggered portion of the AGN population. Thus, 
the maximum possible contribution from major mergers is assumeco The light curve of the 
secularly trigger ed AGN is constrained b y th e minimum to tal % 2 fit to all 12 7 HXLF data 



points . Both the Merloni fc Heinzl (120081 ) and iNetzerl (120091 ) ABHMFs and the lLabita et al. 



( 120091 ) and continuity equation (Eq. \TU^ ABHMF evolutions are considered. It is assumed 
that merger triggered and secularly triggered AGN may have different light curves; thus, 
the light curve parameters for AGN triggered by mergers and AGN triggered by secular 
processes are considered separately. 

We started by attempting to model both the AGN triggered by major mergers and AGN 
triggered by secular processes using the same ABHMF and evolution for both populations 
of AGN. The best fit parameters for these models are summarized by the middle section 
of Table [TJ Despite these models providing decent fits to the observed HXLF data points, 
the models were not in agreement with the other observational constraints, as summarized 
by the middle section of Table [2J However, this exercise did d emonstrate that the z < 0.5 
AGN population is fairly well described using the INetzerl (120091 ) ABHMF with the continuity 
equation evolution (Eq . [TO]) and the z > 1 AGN popu l ation can be described by using the 
Merloni fc Heinzl (boosl ) ABHMF with the lLabita et ah! kood ) evolution. Therefore, a mixed 
ABHMF and evolution model was investigated. 



Since the INetzerl ( 120091 ) ABHMF with continuity equation evoluti on works well at 
l ow re dshift, where secular evolu t ion is expected to dominate, and the iMerloni fc Heinz 
( 120081 ) ABHMF with lLabita et al.l (120091 ) evolution works well at high redshift, where major 
merger triggered AGN are expected to d ominate, a hybrid m odel is calculated. Using the 
Merloni k Heinzl (J2008I ) ABHMF with the lLabita et al.l d2009f ) evolution, the best fit param- 
eters for the merger triggered AGN in this model are to = 2.5 x 10 8 yrs, (5 = 0.7, and rj 



If the maximum contribution from secularly triggered AGN is assumed, an appropriate fit to the HXLF 
can be found (Xred = L7), but this model cannot account for the observed local black hole mass density. 



-14- 



= 2.5. Using the iNetzerl (120091 ) ABHMF with the continuity equation evolution (Eq. ITU]) . 
the AGN triggered by secular mechanisms are best described by the parameters to = 2.5 
x 10 8 yrs, (3 = 0.8, r\ = 0.4, and f sec = 0.01 Gyr -1 . This model does moderately well at 
describing the AGN HXLF, with Xred = 1-4, as shown in Fi gure [21 Also, this mo del is in 
agreement with the local black hole mass density observed bvlShankar et al.l (120091 ) and the 
z ~ 2 black hole mass density observed by lTreister et al.l (|2010bJ), which, as shown in Figure 
[31 is dominated by black hole growth triggered by major mergers. Figure [H shows that this 
model is also in agreement with the observed XRB spectrum. The mixed ABHMF model is 
in decent agreement with the 2-10 keV AGN number counts as shown in Figure [51 Figure [6] 
shows that the mixed ABHMF model is in good ag reement with the Swift/BAT 15-55 keV 
AGN number count observed by lAjello et al.l (120091 ) . Furthermore, at the Swift/BAT survey 



flux limit of 7.6 x 10 erg cm s , we find that merger trig gered AGN accou nt for ~20% 



of 15-55 keV number count, in agreement with the findings of iKoss et al.l ( 120101 ). Thus, this 
mixed ABHMF model provides the best overall fit to the observational constraints. 

The largest discrepancy between the mixed ABHMF model and the observational con- 
straints, is the HXLF at Lx < 10 42 erg s _1 , specifically at z — 0.6. Due to the tendril nature 
of cosmic structure, narrow fields, which tend to provide the high z, low L x AGN detections, 
will observe over dense regions when one of these tendrils passes through the observed field. 
This cosmic vari ance may cause the numb er counts and HXLFs observed in different fields 
to vary slightly (IBrandt fc Hasingerl 120051 ) . Observations of AGN with Lx < 10 42 erg s _1 
and z ~ 0.6 can only be conduc ted in the Chandr a deep fields, which are known to have 
significant field-to- field variation (jCowie et al.ll2002l ). As the slope of the low Lx end of the 
HXLF is still fairly uncertain, cosmic variance is most likely to affect the low Lx end of the 
observed HXLF. 

W hen both merg er triggered and secularly triggered AGN are considered, models which 
use the INetzerl (120091 ) ABHMF are in agreemen t with observations at z < 0.5, but not at 
high redshift. In contrast, models which use the iMerloni fc Heinzl ( 20081 ) ABHMF for both 
merger and secularly triggered AGN over-predict the black hole mass density. In order to 
explain both high and low redshift observations, it is necessary that merger triggered and 
secularly triggered AGN have different ABHMFs which evolve differently with redshift. The 
mass distribution of black holes triggered by secular processes is, therefore, intrinsically 
different than the mass distribu t ion of b lack holes triggered by mergers, in agreement with 
the findings of ISchawinski et al.l (j2010bl ). 
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6. Discussion 

An observationally based AGN triggering rate and a theoretical AGN light curve is 
used to model the AGN population. By considering the HXLF at five different redshifts, the 
black hole mass density, the XRB spectrum, and the AGN number counts in the 2-10 and 
15-55 keV bands, it is determined that neither merger triggered AGN nor AGN triggered 
by secular processes can alone account for the entire AGN population. Two populations are 
necessary to describe the observed evolution and space density of AGN. Furthermore, it is 
found that the ABHMF of merger triggered AGN evolves more rapidly and has a higher 
average black hole mass than the ABHMF of secularly triggered AGN. 



6.1. AGN Light Curve Model 

The AGN light curve model used has three parameters, to, /3, and r]. The r] parameter 
is related to the peak Eddington ratio of the light curve. The AGN population model is least 
sensitive to rj, however, it is found that r\ > 1.0 provides the best fit to the observational 
constraints, suggesting that short periods of super-Eddington accretion are likely in a signif- 
icant fraction of the AGN population. The AGN population model is moderately dependent 
on the time scale parameter, to- If this time scale is too short, the AGN fade too quickly, 
and the space density of low luminosity AGN is severely under-predicted. Thus, t is most 
important for the normalization of the AGN HXLF. The shape of the HXLF is primarily 
controlled by the slope of the light curve, which is controlled by the j3 parameter. The AGN 
population model is most sensitive to /3. Changing /3 by 0.1 can cause the resulting HXLF 
shape to change substantially. The primary goal of this study is to investigate if it is neces- 
sary that AGN be triggered by both mergers and secular processes or if a single mechanism 
can account for the entire AGN population, not to determine the best fit AGN light curve 
parameters. However, the observational constraints used in this study do provide interesting 
restrictions on the AGN light curve parameters. 



Hopkins fc Hernquistl (120091 ) point out that different models of the AGN light curve 
predict different light curve slopes, /3. Therefore, the AGN light curve fit provides information 
about the physical conditions in the vicinity of the accretion disk. For example, self-regulated 
models, where the black hole accretion is feedback limited, predict /3 = 0.3 - 0.8. Meanwhile, 
if AGN are fueled b y mass loss from a nuclear star cluster, /3 = 0.9 - 1.0. Models in which 



the fuel supply of a IShakura fc Sunyaevl ( 119731 ) accretion disk is suddenly cut off, predict /3 



= 0.80 - 0.84. In all of the best fit merger triggered AGN models presented here, /3 = 0.7, 
suggesting that for merger triggered AGN, the accretion is feedback limited. Also, for the 
best fit model presented here, 77, the parameter related to light curve peak Eddington ratio, 
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is much smaller for AGN triggered by secular mechanisms than for AGN triggered by major 
mergers. 



6.2. Two Populations of AGN 

Here, it is assumed that every major merger triggers an AGN and major mergers can 
account for the AGN population at z = 2.3, thus this analysis is an upper limit for the 
contribution of major merger triggered AGN. Still, it is found that at z < 2, a significant 
fraction of all AGN are triggered through secular mechanisms. Figure [7] shows the space and 
luminosity density of AGN with \ogLx > 42, 43, and 44, which are triggered by mergers 
and secular processes as a function of redshift for the mixed ABHMF model. At all three 
luminosity ranges, the luminosity density and space density of AGN is dominated by major 
merger triggered AGN at z > 1.5 and secularly triggered AGN at z < 1.5. By varying 
model parameters, we find the minimum redshift at whic h major mergers c a n dom inate 
AGN trigger is z ~ 1. This finding is in agreement with iGeorgakakis et all (120091 ). who 
found that ~30% of the AGN space density and ~25% of the AGN luminosity density at 
z ~ 1 is due to AGN hosted by disk dominated hosts, implying that secular processes are 
responsible for at least one quarter of the AGN luminosity density at z ~ 1. It is found that 
merger triggered AGN account f or ^20% of the sp ace density of log L x > 43 AGN at z < 
0.05, consistent with findings of iKoss et al.l ( 12010I ). Thus, AGN hosted by disk dominated 
galaxies are a significant fraction of the AGN population by z ~ 1, indicating that secular 
evolution is an important mode of galaxy evolution at this redshift. 



Similarly, iDraper fc Ballantynd (j2011bl ) found that for AGN at z < 1 the host galaxies 
of obscured and unobscured AGN are, on average, the same, suggesting that at z < 1 AGN 
and galaxy evolution is controlled by secular processes. However, at z > 1, the host galaxies 
of unobscured AGN are intrinsically less dusty than the host galaxies of obscured AGN at 
similar redshifts, s uggesting a fundamental cha nge in the mechanisms which control AGN 
activity at z ~ 1 ( IDraper fc Ballantynd 1201 lbl ). Indeed, Figure [2] shows that, at z < 1, 
AGN triggered through secular processes dominate even the high Lx end of the HXLF, 
while fading major merger triggered AGN dominate the low Lx end of the HXLF. These 
findings are also in agreement with the conclusions of a variety of recent observational studies 
which found that secular proc esses are an important form, and possibly the dominate form, 



of ga l axy evolution at z < 2 (IGeorgakakis et al. 



2011; Schawinski et al. 2011: Williams et al 



2009: Allevato et al. 2011; Cisternas et al. 



l201ll ). Indeed, according to the HXLF, we find 



that AGN triggered by secular processes can a ccoun t for the entire AGN population with 
logLx > 43 at z < 0.5. iHopkins &: Hernquistl ( 120061 ) suggest that secularly triggered AGN 
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are not cosmo logically important; however, the findings of this study illustrate that AGN 
triggered by secular mechanisms are necessary to describe the AGN population and dominate 
the space density of AGN with logLx > 43 at z < 1.5. 



The findings of ISchawinski et al.l ( 12 01 Obi ) suggest that this fundamental change in AGN 
activity at z ~ 1.5 is due to cosmic downsizing. By studying AGN black hole masses as 
a functi on of host galaxy morpho logy for AGN observed by the Sloan Digital Sky Survey 
(SDSS), ISchawinski et al.l (j2010bl ) find that it is preferentially the least massive black holes 
in less massive early type galaxies (stellar mass M* ~ 10 10 M Q ) which are currently active. 
In contrast, the black holes which are currently active in late type galaxies are preferentially 
the most massive black holes hosted by more massive late type host galaxies (M* ~ 10 11 
M ). A stellar disk generally indicates that a galaxy has not undergone a major merger 
or that the last major merger e xperienced by the galaxy was long enough ago that a disk 
had time to reform, > 500 Myr (IHota et al.ll20f 11 ). We can therefore assume that, for AGN 
hosted by late type galaxies , the cu rrent episode of AGN activity was triggered by secular 



processes. ISchawinski et al.1 (j2010bl ) suggest that the early type galaxy hosts in their sample 
may be smaller versions of the mode of galaxy evolution experience d by massive ellip t icals at 
high redshift. Thus, the early type host galaxies in the sample of ISchawinski et al.l (J2010bJ) 
may host downsized versions of major merger triggered AGN. Thus, it appears that the 
dominance of major mergers in galaxy evolution began to decline by z ~ 2 and is continuing 
to decline, in agreement with the theory of cosmic downsizing. Meanwhile, the importance 
of secular evolution increased as the importance of major mergers decreased, until secular 
evolution became the dominate form of galaxy evolution at z ~ 1.5. Since we are still 
early in the era of secularly dominated galaxy evolution, it is the more massive systems 
which are currently undergoing secularly triggered AGN activity. Cosmic downsizing would 
therefore predict that as the era of secular dominance continues, AGN triggered through 
secular processes will become more common in lower mass systems. 



6.3. ABHMF and Evolution 



It is found that AGN triggered by major mergers at z < 4 have, on average, more massive 
central black holes than AGN triggered by secular processes, at least at high redshift. Locally, 
it appea rs that merger triggered AGN and secularly triggered AGN have similar black holes 
masses (ISchawinski et al.ll2010bl ). necessitating that the ABHMF of merger triggered AGN 
evolves differently than the ABHMF of secularly triggered AGN. The ABHMF of AGN 
triggered by secular mechanisms appears to evolve in a manner consistent with the continuity 
equation (Eq. HD|) . which describes the evolution of the combined active and quiescent black 
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hole mass function (e.g., ISmall fc Blandfordlll992l ; iMerloni &: Heinzl I2008J ) . However, if the 
ABHMF of merger triggered AGN evolves following the continuity equation, then at z > 1 
the average active black hole mass is t oo small to account for the high Lx end of the HXLF. 
In contrast, if the lLabita et all (120091 ) ABHMF evolution is used to evolve the ABHMF of 
secularly triggered AGN, then the low Lx end of the z < 1 HXLF is significantly under- 
predicted. Thus, it appears that not only are both secular processes and major mergers 
important mechanisms for triggering the z < 2 AGN population, but these two mechanisms 
trigger different populations of black holes and these two populations have different redshift 
evolution. 

It is not surprising that the ABHMF of merger triggered AGN evolves differently from 
the secularly triggered ABHMF and the quiescent black h ole mass function. T he major 
merger rate evolves strongly with redshift and galaxy mass (IHopkins et al.ll2010al ). Thus, it 
is expected that the ABHMF of AGN tr iggered by major mergers would also evolve strongly 
with redshift (e.g., iTreister et al.ll201ll ). In contrast, the rate of AGN being triggered by 
secular evolution, f sec , is assumed to be constant with redshift. Therefore, it is expected 
that the ABHMF for AGN triggered by secular processes will be less redshift dependent than 
the ABHMF for merger triggered AGN. Also, as the majority of supermassive black holes 
at any redshift are quiescent, the combined quiescent and active black hole mass function 
evolves more slowly with redshift than the merger triggered ABHMF. Thus, it is expected 
that the merger triggered ABHMF will evolve more strongly with redshift than the ABHMF 
of secularly triggered AGN or of the combined active and quiescent black hole mass function. 



6.4. The Connection Between Compton Thick AGN and Mergers 



Observational and theoretical evidence suggests that Compton thick AGN are AGN 
which were recently triggered by a mechanism which has caused a large a mount of gas 
and dust to be 



Page et al. 



2004 



unneled into the nuclear regions of the host galaxy (e.g., Fabian 



Fabian et alJ 12008 . 120091 : Draper fc BallantyndbOlOt ITreister et al.ll2010al : 



1999 



Draper fc Ballantynd l2011bl ). Due to the large amount of gas and dust required for the 
rapid AGN fueling and high column density necessary for an AGN to be Compton thick, 



the mo st likely triggering process is a wet major merger (ISanders et al.lll988t IHopkins et al. 



2006aj), though I Weinzirl et al] (120111 ) argue that Compton thick levels of obscuration are also 
possible in instances of clumpy cold-flow accretion onto the host galaxy. Here, we set for 
= 0.0 for the secularly triggered AGN and calculate fcr for the merger triggered AGN as 
described in Section 13.11 If it is assumed that Compton thickness is an evolutionary phase 
which only merger triggered AGN experience, the Xred °f the mixed ABHMF increases by 
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0.1 and the model is still in agreement with the observed black hole space density, XRB 
spectrum, and the 2-10 keV and 15-55 keV number counts. Therefore, the hypothesis that 
Compton thick AGN are an evolutionary stage of merger triggered AGN is fully consistent 
with the model presented here. 



7. Summary 

By combining an observationally motivated AGN triggering rate and a theoretically 
motivated AGN light curve, an AGN population synthesis model is developed which can track 
the evolution of AGN triggered by mergers and secular processes. The light curve parameters 
are constrained by comparing model predictions against the observed AGN HXLF. The 
observed space density of supermassive black holes, AGN number counts, and XRB spectrum 
are also considered. The results of this study are summarized below. 

• Neither AGN triggered by major mergers nor AGN triggered by secular mechanisms 
alone are capable of accounting for the entire observed AGN population at z < 2. This 
finding is independent of the SMF used or the potential AGN host galaxy minimum mass. 

• The dominant AGN triggering mechanism switches from major mergers at z > 1.5 
to secular mechanisms by z ~ 1. This is a natural consequence of both the major merger 
rate and fraction of gas rich galaxies increasing with redshift. However, both major mergers 
and secular processes are important triggering mechanisms for AGN over the entire redshift 
range considered, especially at moderate luminosities. . 

• The black hole mass density is dominated by black hole growth triggered by major 
mergers over the entire redshift range considered. 

• The space density and luminosity density of AGN with logLx > 42, 43, and 44 is 
dominated by AGN triggered by major mergers at z > 1.5 and by AGN triggered by secular 
processes at z < 1.5. 

• The XRB is dominated by emission from AGN triggered by secular mechanisms. 

• The ABHMF of AGN triggered by major mergers is different from the ABHMF of 
AGN triggered by secular mechanisms. At high redshift, the average black hole triggered by 
major mergers is more massive than the average black hole triggered by secular processes, 
which is consistent with cosmic downsizing. 

• The evolution of the ABHMF of AGN triggered by major mergers is also different 
from the evolution of the ABHMF of AGN triggered through secular mechanisms, with 



-20- 



the ABHMF of major merger triggered AGN evolving more rapidly with reds hift than the 



ABHM F of AGN triggered by secular processes. The observed evolution of lLabita et al. 



( 120091 ) appears to describe the evolution of the merger triggered AGN ABHMF well. In con- 
trast, the continuity equation (Eq. [TO]) describes the evolution of the ABHMF for secularly 
triggered AGN. This finding is consistent with the strong redshift evolution of the major 
merger rate. 

• The findings of this study are not changed if Compton thick AGN are an evolutionary 
stage only experienced by major merger triggered AGN. 

This work was supported by NSF award AST 1008067. The authors thank the referee 
for helpful comments which improved this paper. 
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Table 1. Summary of derived model free parameters. 



Model Merger Triggered AGN Secularly Triggered AGN 

*o (yr) P V to (yr) P V fsec (Gyr -1 ) 



One Population Models 



Mergers only 
Secular only 



2.5 x 10 8 0.7 2.5 



2.5 x 10 8 0.8 0.4 



0.0 
0.02 



Two Population Models 



N09 w/ Eq. EH 2.5 x 10 8 0.7 0.4 

N09 w/ L09 1 x 10 8 0.5 1.0 

MH08 w/ Eq. \W\ 2.5 x 10 8 0.8 2.5 

MH08 w/ L09 2.5 x 10 8 0.7 2.5 



5 x 10 8 0.5 2.5 

2.5 x 10 8 0.6 2.5 

5 x 10 8 0.7 2.5 

5 x 10 8 0.8 2.5 



0.005 
0.005 
0.01 
0.005 



Two Population Mixed ABHMF Model 



Mixed ABHMF 



2.5 x 10 8 0.7 2.5 I 2.5 x 10 8 0.8 2.5 



0.01 



Note. 



MH08 refers to the 



ABHMF, and L09 refers to the Labita et al. 



Merloni fc HeinzJ (hoosh ABHMF, N09 refers to the iNetzeil (hoOfll ) 



(J2009|) ABHMF evolution. The mixed ABHMF model 
uses the MH08 ABHMF with L09 evolution for the merger triggered AGN and the N09 ABHMF 
with the continuity equation (Eq. [T0|) evolution for the secularly triggered AGN. 
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Table 2. Summary of model fits to the observational constraints. 



Model 


d -°- f - a X 2 red b P. c XRB d 


2-10 keV e 15-55 keV f 




One Population Models 


Mergers only 
Secular only 


123 2.4 y x 
122 2.1 x x 


X X 
X X 


Two Population Models 



N09 w/ Eq. [10] 119 1.5 x 

N09 w/ L09 119 1.4 x 

MH08 w/ Eq. M 119 1.7 x 

MH08 w/ L09 119 1.5 j 



X 


V 


V 


X 


V 


V 


X 


X 


X 


X 


X 


X 



Two Population Mixed ABHMF Model 



Mixed ABHMFS 119 1.4 



V V V V 



a The number of degrees of freedom for the x 2 test reported in column Xred- 

b The quoted x 2 e d re f ers to the 127 HXLF data points plotted in Figure [TJ 

C A x in this column denotes that, when considering the two black hole space density data points, 
X 2 > 4.0 for the designated model. 

d A x in this column denotes that this model either over or under-predicts the XRB by more 
than 25% in the 10-50 keV band. 

C A x in this column denotes that, this model either over or under-predicts the the 2-10 keV 
observed number counts shown in Figure [5] by more than 25% at all fluxes. 

A x in this column denotes that, this model either over or under-predicts the 15-55 keV observed 
logA r -log5 shown in Figure O by more than 25% at all fluxes. 

g The mixed ABHMF model uses the MH08 ABHMF with L09 evolution for the merger triggered 
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AGN and the N09 ABHMF with the continuity equation (Eq. [T0|) evolution for the secularly 
triggered AGN. 
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Fig. 1. — Best fit one population models compared against the observed HXLF. The dot- 
dashed red lines show the HXLF from the best fit major merger trigger only model and the 
dashed blue lines show the HXLF from the best fit secular processes trigger only model. The 
data points show measurements of the HXLF by Ueda et al. (2003; filled circles), La Franca 
et al. (2005; stars), Silverman et al. (2008; triangles), Aird et al. (2010; squares), and Ueda 
et al. (2011; open circles). 
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Fig. 2. — The HXLF predicted by the mixed ABHMF model, which has Xred = 1-4- The 
solid black line is the total AGN HXLF. The dot-dashed red lines show the contributions 
from AGN triggered by mergers and the dashed blue lines show the contribution from AGN 
triggered by secular mechanisms. The data points are the same as in Figure [TJ 
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Fig. 3. — The black hole mass density with respect to redshift, predicted by the mixed 
ABHMF model. The line styles are the sa me as in Figur e [2] Th e z ~ shaded region is the 
local black hole mass density observed b y IShankar et al.l (120091) and the z ~ 2 data point is 
the black hole mass density observed by iTreister et al.l ( 12 01 Obi ) . 
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Fig. 4. — The XRB spectrum predicted by the mixed ABHMF model. The line styles are 
the same as in Figure EJ Colored regions and data points show measurements from various 
instruments; blue: ASCA GIS (Kushino et al. 2002); magenta: Rossi X-ray Timing Explorer 
(RXTE; Revnivtsev et al. 2003); green: XMM- Newton (Lumb et al. 2002); red: BeppoSAX 
(Vecchi et al. 1999); yellow: ASCA SIS (Gendreau et al. 1995); cyan: XMM- Newton (De 
Luca & Molendi 2004); grey data: HEAO-1 (Gruber et al. 1999); blue data: INTEGRAL 
(Churazov et al. 2007); red data: Swift/B AT (Ajello et al. 2008); black data: Swift/XRT 
(Moretti et al.2009); green data: INTEGRAL (Tiirler et al. 2010). 
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Fig. 5.— The 2-10 keV AGN JV(> S) predicted by the mixed ABHMF model. The line styles 
are the same as in Figur e [21 The plotted d ata show observed t he number counts from various 



surveys: C-C OSMOS (Elvis et al. 



and XMM-Newton Serendipitous (IMateos et all 120081 ). 



20091 ). XMM-COSMOS (ICappellutil l2009h. Swif t/XRT 



Serendipitous (IPuccetti et al.1120111). XMM-Newton Lockman Hole ( IBrunner et al.l 120081 ). 
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Fig. 6.— The 15-55 keV AGN JV(> S) predicted by the mixed ABHMF model. The line 
styles are the same as in Figure P The data point s show the 15-55 keV AGN \ogN-\ogS 
relation observed by Swift/B AT JAiello et allbopgh . 
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Fig. 7. — The space and luminosity density of AGN as a function of redshift predicted by 
the mixed ABHMF model. The top row shows the space density of AGN for logLx > 42, 

43, and 44. The bottom row shows the luminosity density of AGN for logLx > 42, 43, and 

44. The line styles are the same as in Figure [H 



